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We present two novel approaches to establish the local 
density of states as an order parameter field for the Anderson 
transition problem. We first demonstrate for 2D quantum 
Hall systems the validity of conformal scaling relations which 
are characteristic of order parameter fields. Second we show 
the equivalence between the critical statistics of eigenvectors 
of the Hamiltonian and of the transfer matrix, respectively. 
Based on this equivalence we obtain the order parameter ex- 
ponent ao ~ 3.4 for 3D quantum Hall systems. 



The absence of diffusion in coherent disordered elec- 
tron systems is known as Anderson localization [Q. 
In dimensions d > 2 a disorder induced localization- 
delocalization (LD) transition occurs quite generally at 
some value of the Fermi energy || . In d — 2 all states are 
localized unless a certain amount of spin-orbit scattering 
U or a strong magnetic field is present [||||. The LD 
transitions of independent (spinless) 2D electrons sub- 
ject to a strong perpendicular magnetic field are located 
at the Landau energies. These transitions are generally 
believed to be responsible for the integer quantum Hall 
effect §,§. 

In general, LD transitions are characterized by the 
critical exponent v of the localization length M and 
by the multifractal /(a) spectrum of the local ampli- 
tudes of critical eigenstates The /(a) spectrum 
describes the statistics and scaling behavior of the local 
density of states (LDOS). Although the average density 
of states does not reflect the LD transition, the typical 
value (i.e. geometrically averaged) of the LDOS does: it 
vanishes with an exponent /3 typ = (ao — d)v on approach- 
ing the LD transition point, where ao > d is the maxi- 
mum position of /(a). It is thus tempting to interpret 
the LDOS as an order parameter field of the LD problem 
H (see also |n]|ll|] ) ■ Our aim is to support this interpre- 
tation by establishing two characteristic features of order 
parameter fields for the LDOS: First, scaling exponents of 
the order parameter field are related to the critical expo- 
nents of the corresponding spatial correlation functions. 
These correlation functions show conformal invariance. 
Second, the scaling exponents are universal in the sense 
of (one-parameter) scaling theory, i.e. any local quantity 
containing contributions from the relevant scaling field 
shows asymptotically the same spectrum of scaling ex- 
ponents. 



In this article we demonstrate that the /(a) spectrum 
of critical eigenstates is related to correlation functions in 
different geometries by conformal invariance. We derive 
the conformal mapping relations appropriate for a multi- 
fractal situation and check them numerically in 2D quan- 
tum Hall systems (QHS). Furthermore, we have calcu- 
lated numerically /(a) for the local components of trans- 
fer matrix eigenvectors in 2D QHS and show that it co- 
incides with /(a) of the Hamiltonian eigenstates. Thus, 
these two local quantities share the same spectrum of 
scaling exponents although their microscopic construc- 
tion is quite different. Our findings support the identifi- 
cation of the LDOS as an order parameter field for the 
LD transition. From /(a) calculations for the transfer 
matrix eigenvectors in a 3D QHS we obtain the charac- 
teristic order parameter exponent ao ~ 3.4. 

The multifractal analysis of the statistics of critical 
eigenstates ip (see e.g. H) usually starts from considering 
the box-probability, P{L b ) ■— J hox bA| 2 j 01 some box with 
linear size L b . The box- probability is normalized to the 
total volume L d , P(L) = 1. At the LD transition the 
corresponding distribution function 7r(P; L b /L) gives rise 
to power law scaling for the moments 



(P q {L b )) L cx {L b /L) d+T ^ 



(1) 



where r(q) is non-linear. The distribution function can 
be described in terms of a single-humped, positive, func- 
tion f(a): 



ir(P; L b /L) dP oc (Z 6 /£) d ~ /(a) da . 



(2) 



where a := \n P / hi(L b / L) . The multifractal spectrum 
/(a) is the Legendre transform of r(q). Thus, the 
statistics of critical eigenstates is encoded in /(a) or 
equivalently in r(q). The statistics of critical eigen- 
states transforms to the statistics of the LDOS p(r) = 
A _1 |?/;(r)| 2 at criticality. Here A denotes the mean level 
spacing at the corresponding energy. It is a non-critical 
quantity, A cx L~ d . Consequently, the typical value, 
pty P := exp(ln ( o) oc L d ~ a ° where ao is the maximum 
position of /(a). Scaling L with the localization length 
shows that pt yp vanishes on approaching the transition 
point with characteristic exponent /3 typ = (ao — d)v. This 
observation motivates to consider the LDOS as an order 
parameter field of the LD transition || . 

To make contact between /(a) and correlation expo- 
nents we first make the general scaling ansatz for the 
q-dependent spatial correlation in square-like geometries 
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(p q (v)p q (v')) L oc a(r/L) ■ r" 5 ^ , r = 
Here z(q) fulfills the scaling relation 

z(q)= 2(1 -q)d + 2r(q). 



(3) 



(4) 



In ordinary critical phenomena the coefficient a(r/L) 
will approach a constant for r <C L — > oo. Then, stan- 



dard conformal mapping arguments in 2D 13 say that 
the corresponding correlator in a long strip of width Lt 
with periodic boundary conditions in transverse direction 
decays like 



exp[-2u/^(Lr)], 



(5) 



where u is the distance in longitudinal direction and 
^(Lt) define generalized (g-dependent) localization 
lengths in quasi-lD. The localization lengths are then 
related to r(q) in the following way § (see also @) 

e q] (L T yL T = 2[7rz(q)]- 1 . (6) 

Equation (^J) establishes a relation between generalized 
localization lengths and f(a). 

In the following we will show that this relation remains 
valid, although the coefficient a(r/L) in Eq. (||) cannot 
be treated as a constant. In a multifractal situation the 
coefficient a(r/L) in Eq. (|^) remains sensitive to the ac- 
tual system size L even for L — > oo (at the transition 
point). It has been shown that a(r/L) oc (r/L) V2 ^ q > with 
2/2(9) = d + r(2q) — 2qd [p||il|. Thus, standard conformal 
mapping arguments do not apply. This fact has been 
stressed recently by Zirnbauer |^5|. We recall that the 
main idea behind conformal mapping arguments is an 
extension of a homogeneity law for correlation functions 
with respect to rescaling. Such a law exists also in our 
case, since any rescaling of all length scales in the corre- 
lator of Eq. (0) by the same scaling factor s leads to 



(p*( S r)^( S r')> 



sL 



(p'(r)p 9 (r'))i 



(7) 



We make now the plausible assumption that this law can 
be extended to conformal mappings which are local scale 
transformations: conformal mappings transform the cor- 
relator of a geometry f2 to the corresponding correlator 
in geometry f2. For large but finite 2D systems we are 
led to 



(p q (w( Zl ))p q (w(z 2 ))) {2 
(p q (z 1 )p q (z 2 )) n 



^ \w'(z 2 )\ 



(8) 



where w(z) is any holomorphic function of complex co- 
ordinate z and w'(z) denotes the derivative. 

By choosing w(z) — (Lt /2tt) In z which maps the plane 
onto a strip straightforward calculations (cf. |ll|) show 
that the correlator in the strip is given by expression rtq) 
with an additional prefactor, exp[(2iri/2(q) / Lt)(u — L s )\, 
where L s is the strip length. Note, that the important 
relation Eq. (^|) between ^ and z(q) remains valid. 



To calculate the generalized localization lengths £^1 in 
strip- like (quasi- ID) systems (L s being the length and 
Lt -C L s being the width in the remaining d — 1 direc- 
tions) one can start from a Hamiltonian H. The Green's 
function G(r,r';E) := (r'| (E - H + £0+) _1 |r) can be 
calculated for r, r' situated at opposite ends of a strip by 
a recursive method |l6| . The resulting Green's function 
G(L S ;E) yields (cf. expression (0)) 



-(2L s )- 1 ln(\G(L s )\ 2q ) =: \® . 



(9) 



By Eq. (||), A^] is proportional to z(q) which becomes 
negative for q > 1. This corresponds to those rare events 
where the Green's function increases with L s . To detect a 
reasonable amount of these rare events one has to restrict 
the numerical calculation to strip-lengths L s of a few 
typical localization lengths. 

Alternatively, one can construct a transfer matrix T 
which relates either the values of the wavefunction at op- 
posite ends of the strip or relates the corresponding scat- 
tering states Jl7| , |T8|] . For a transverse width Lt — Ml, I 
being a microscopic scale of the problem, the dimension 
of the transfer matrix T is proportional to M. Lyapunov 
exponents, pt, can be defined from the eigenvalues of TT^ 
corresponding to eigenvectors u^, 



(10) 



The smallest positive Lyapunov exponent, denoted by 
p, determines the generalized inverse localization lengths 

ira, 



A^ = -(2L s )- 1 \n(e- qli ) 



(11) 



We have calculated A^ for two models of QHS. The 
first model is the random Landau model (RLM) describ- 
ing the Hamiltonian of 2D disordered electrons restricted 
to the Hilbert space of one Landau band. The disorder 
is caused by a white-noise potential. The characteristic 
microscopic scale is the magnetic length l},. A detailed de- 
scription of this model is given e.g. in Ref. j|. We applied 
the recursive Green's function method (RLM/GF) and a 
transfer matrix method (RLM/TM). The second model 
is the network model (NWM) of Chalker and Codding- 
ton 1 18 1 (with a characteristic microscopic scale I) which 
is formulated in terms of a transfer matrix. 

In Fig. 1. we show the corresponding X^Lt curves and 
its first derivative for the NWM in the regime \q\ < 1. 
The conformal mapping relation, Eq. (|^), allows for the 
calculation of f(a) and, especially, of the most interesting 
scaling exponent ckq. For comparison with related work 
|p(if we also give the values for r(2). We find within 
the RLM by the GF method a = 2.28 ± 0.03 for strip 
widths of lOO^s and by the TM method a = 2.3 ± 0.03 
and t(2) = 1.66 ± 0.03 for strip widths of 50/^. For the 
NWM we obtain a = 2.27 ± 0.01 and r(2) = 1.64 ± 
0.03 for strip widths of 64/. The results agree for the 
different models and agree with the values reported in 
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the literature |20| . We emphasize that our results confirm 
the conformal mapping relation, Eq. ([|). 

In three dimensions only a few conformal mappings 
exist. Especially, for three dimensional strips with pe- 
riodic boundary conditions we have no appropriate con- 
formal mapping at our disposal. We saw that conformal 
invariance is sufficient to derive relations between f(a) 
and \w , however it may not be necessary. Therefore, 
we found it worthwhile to study if a similar relation as 
Eq. (||) could be valid for 3D strip-like systems. We thus 
make the ansatz 

L T \^{L T ) = C{r{q)- 3(<? - 1)) , (12) 

where C is some unknown constant. Since r(q) has to 
vanish at q = 1 and is smaller than d(q— 1) for q > 1 this 
ansatz raises the question if the generalized localization 
lengths A^l vanish for q = 1 and become negative for 
q > 1, for the 3D problem too. To check for such behavior 
we took a 3D version of the network model, introduced 
recently by Chalker and Dohmen |2l}| . The 3D network 
model consists of coupled layers of 2D networks and the 
corresponding transfer matrix is described in detail in 
Ref. pl| . Our numerical findings at LD transition points 
of the model are shown in Fig. 2. They indicate that 
A^l is compatible with our ansatz. Note that, so far, we 
have established a relation between f(a) and A^ only at 
criticality in 2D. Our 3D results suggest that a similar 
relation might be valid in 3D at criticality. To us, it 
seems rather unlikely that similar relations hold true far 
off criticality. 

Our second approach to establish the LDOS as an or- 
der parameter for the LD transition is based on the as- 
sumption that any local quantity in a given disordered 
system can be viewed as being decomposed in terms of 
local scaling fields in the sense of critical phenomena the- 
ory (e.g. p2| ). If one-parameter scaling holds true then, 
on approaching criticality, only the relevant scaling field 
has to be considered giving rise to a universal distribu- 
tion of local amplitudes. Any contribution of irrelevant 
scaling fields would lead to deviations between the cor- 
responding f(a) spectra, but eventually die out in the 
thermodynamic limit. 

The components of the eigenvectors corresponding 
to the transfer matrix, Eq. (|l(j) , in the network model |lq | 
have local meaning (the corresponding scattering chan- 
nels are defined in real space) and can be subject to a 
conventional multifractal analysis. From this one obtains 
a one-dimensional multifractal spectrum. Multiplication 
with d yields the associated d-dimensional spectrum. 

When considering the components of transfer matrix 
eigenvectors which belong to large Lyapunov exponents 
(small localization lengths) deviations from a universal 
scaling behavior can clearly be observed in the 2D net- 
work model. For the smaller Lyapunov exponents, how- 
ever, we obtain almost identical results (see Fig. 3). Con- 
centrating on the eigenvector u corresponding to the 



smallest Lyapunov exponent, we find ao = 2.3 ± 0.02 
and t(2) = 1.59 ± 0.05. For the 2D network model we 
have thus evidence that eigenvector components of the 
transfer matrix give rise to the same f(a) spectrum as 
the eigenstates of the Hamiltonian. 

We also calculated /(a) of transfer matrix eigenvec- 
tors for the 3D network model at LD transition points. 
It is worth mentioning that the 3D network model is 
anisotropic and for most values of the parameter that 
fixes the coupling between adjacent layers the typical 
quasi- Id localization lengths are smaller than the system 
width Lt in one of the transverse directions, but the ba- 
sic condition for the applicability of the boxing method is 
that the localization length has to be larger than Lt- Re- 
stricting to those values of the coupling where this condi- 
tion was fulfilled we find from the statistics of the eigen- 
vector components corresponding to the smallest Lya- 
punov exponent: a = 3.4 ±0.2 and r(2) = 2.3 ±0.2 j23|. 
We mention that these values of ao and t(2) fulfill the 
inequality 2/t(2) < v < (ao — d)^ 1 proposed by Janssen 
|S| relying on one-parameter scaling arguments (v « 1.35 

In summary, we have proposed a conformal scaling re- 
lation, Eq. (](]), between the f(a) spectrum of the local 
density of states at criticality and the corresponding gen- 
eralized Lyapunov exponents in strip-like geometries. We 
have checked its validity for 2D quantum Hall systems 
by numerical calculations. We have presented numer- 
ical results which support the expectation that similar 
relations are valid in 3D systems at criticality. Further- 
more, our numerical calculations show that eigenvector 
components of the transfer matrix give rise to the same 
f(a) as that of the Hamiltonian eigenstates. We like to 
stress the following aspects of our results: First, they 
show that the local density of states has several features 
in common with order parameter fields in ordinary crit- 
ical phenomena. Second, the very existence of relations 
between the critical statistics of eigenstates (or LDOS) 
of the Hamiltonian, the generalized Lyapunov exponents, 
and the critical statistics of eigenvectors of the transfer- 
matrix is new and far from obvious. Third, such relations 
provide alternative ways of computing the f(a) spectrum 
of critical eigenstates. As an application of this we calcu- 
lated the critical order parameter exponent for 3D quan- 
tum Hall systems. 
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Figure 1: Generalized inverse localization lengths X^Lt 
(A) and first derivative with respect to q (B) for different 
strip widths Lt — Ml in the 2D network model (I is the 
corresponding microscopic length scale). 

Figure 2: Generalized inverse localization lengths A^Lt 
(A) and first derivative with respect to q (B) for different 
cross-sections (Lt) 2 = (Ml) 2 in the 3D network model 
(I is the corresponding microscopic length scale). 

Figure 3: The exponent ao describing the statistics of 
eigenvector components corresponding to the z-th Lya- 
punov exponent in the 2D network model. 
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